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Abstract. We argue that the interaction of stellar wind with the surface of a circumstellar accretion (or proto- 
planetary) disk can result in the acceleration of relativistic electrons in an external layer of the disk, and produce 
synchrotron radiation. Conservative estimates give a total synchrotron luminosity Ls ~ W~^Lq for a central star 
with AI — 10"'' Mq yr~^, comparable with the value observed around the TW object in the W3(OH) region. 
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1. Introduction 

Recent investigations have shown that H2O masers in the 
vicinity of young stellar objects can be associated both 
with jets and circumstellar disks (see Torrelles et al. 120021 
for the most recent review) . It is suggested that systems in 
which H2O masers trace disks are less evolved than those 
in which masers trace outflows. 

The detection of synchrotron emission a.t v = 1.6—14.7 
GHz from a very young massive stellar object associated 
with H2O masers in the W3(0H) region (Turner & Welch 
1984, hereafter TW object) demonstrated the presence of 
relativistic electrons in the region nearly coincident with 
the location of the H2O masers (Reid et al. 11995(1 . These 
authors interpreted the synchrotron emission as related 
to acceleration of the relativistic electrons in shock waves 
associated with a powerful jet with E-W orientation em- 
anating from the young stellar object (Reid et al. 119951 
Wilner, Reid, & Menten 1999). Association with the jet 
is based on the measurements of proper motions of the 
H2O masers which are compatible with the hypothesis of 
a bipolar outflow moving in the E-W direction (Alcolea 
et al. 11993(1 . We argue that interpretation of proper mo- 
tions pattern of the H2O maser by Alcolea et al. fl993') is 
not unique (a similar situation is described in Fiebig et al. 
11996(1 and does not rule out the hypothesis of formation of 
the H2O masers on the surface of the circumstellar disk or 
in a molecular ring seen edge-on. Evidence pro and contra 
the disk hypothesis are presented throughout paper. 
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Interstellar H2O masers are usually located in regions 
influenced by strong MHD shocks emanating from young 
stellar objects. These shocks affect the boundaries of the 
jets, outflows, the surface of circumstellar disks and the 
parts of the ambient molecular core that are closest to the 
young stellar object. So one might, in principle, expect 
that the acceleration of relativistic electrons generating 
synchrotron emission can arise not only at the jet bound- 
aries and can be a common phenomenon for the vicinities 
of H2O masers. The question therefore arises whether syn- 
chrotron emission at locations close to H2O masers does 
necessarily originate in jets or if an alternative explana- 
tion connected with high energy processes in circumstellar 
disks is possible. 



In this paper we propose estimates which show that 
the interaction of stellar wind from a central star with an 
accretion disk can support favorable conditions for the ac- 
celeration of electrons and the generation of synchrotron 
radiation from relatively narrow boundary layers on the 
disk surface. We show that this interaction can produce 
a total synchrotron luminosity in the TW object com- 
parable to what is observed. In Section 2 we present sim- 
ple estimates of the acceleration mechanism, and calculate 
the resulting synchrotron emission from the disk surface, 
in Section 3 we describe the spatial distribution of the 
synchrotron radiation, in Section 4 we provide arguments 
in favor of a circumstellar disk around the TW object, 
Section 5 summarizes the results. 
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2. Energy requirements 

Strong winds from early-type stars generate through hy- 
drodynamical instabihties numerous shocks with Mach 
numbers up to 20 (Carlberg I198UL Owocki & Rybicki 
I1984|l . The density of shocks increases near the disk sur- 
face due to the Kelvin-Helmholtz (KH) instabihty of the 
shear flow between the wind and the surface, and subse- 
quent interaction of the wind with the vortices and coher- 
ent structures generated by the instabihty. On nonHnear 
stages these structures have sizes and spacing of the order 
of the thickness of the mixing layer (Roshko 119761 Head 
& BandvoDadhvav [T981|l : 5 ^ y/^rjrjv^. Assume that the 
energy density of the shocks in a layer of thickness AH{r) 
is a fraction f3 of the wind energy density: £s = ppv'^/2, 
where 



M 



(1) 



Woo is the terminal velocity of the wind, M , the mass loss 
rate. The total energy of shocks in the boundary layers 
from both sides of the disk is 



f3Mvo 



R 



AH{r) 



dr, 



(2) 



where is the stellar radius, R, the disk radius. Assume 
that AH{r) = hr, ft.< 1. Then 



= 10^^(3hM^eVsR3 erg. 



(3) 



where M_e = M/{10-^Mq yr^i), = «oo/(10^ cm s^^), 
i?3 — i?/(10'^ AU) (note that the radius of the synchrotron 
structure observed by Reid et. al. 119951 is of the order of 
103 AU). 

Charged particles will be accelerated by the shocks 
through the first order Fermi mechanism. Assume that 
the mechanical energy of shock waves is transformed into 
cosmic ray energy until the nonlinear back reaction of ac- 
celerated particles leads to significant dissipation of the 
large-scale turbulence (Zank, Axford, & McKenzie 1990). 
Under these conditions one can expect that in the steady 
state the energy accumulated in non-thermal ions will 
be comparable to the energy stored in shocks, while the 
fraction of energy contained in relativistic electrons is 
(rrie/mp)^''"^^/^, q being the exponent of a power law spec- 
trum of relativistic electrons (see below). The maximum 
energy of non-thermal ions is found from the condition 
(Bykov & Fleishman 1992) 



3uls 
vA 



1, 



(4) 



where u is the shock speed, Is is the distance between the 
shocks, V is the particle velocity, A is the transport length 
(Toptygin 1985, Bykov & Fleishman 1992) 
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where — pc/eB, ^ is the power-law exponent of the 
fluctuating magnetic field spectrum d{b^)/dk ^ k^^ at 
kls ^ 1 (Bykov & Fleishman 1992). The mean distance 
between the shocks in strong radiation-driven winds is of 
the order Is - lO-^R^{r/R^)^ (Lucy 1982, White 1985), 
and can reach Is ^ 0.1 — 10 AU at r ~ 1 — 10 AU for an 
O- or early B-type star. The actual value can be smaller 
if multiple reflections of the shocks on the vortices and 
coherent structures in the boundary layer of the disk are 
taken into account. As a conservative estimate one can 
take - 1 AU. This results in ^ 300 GeV for ^ = 
1.5, and e,„ ~ 10 GeV for ^ = 1.7. We will assume that 
the energy spectrum of non-thermal electrons extends also 
to these values, so that the characteristic frequency I'b — 
2TT'yeB /rrieC^ , where 7 — e/meC^, can be as high as 20 
GHz. 

With these assumptions and for a power-law spectrum 
of relativistic electrons 

N(e)de — ne^'^de, e > e,nin, 9 > 2, (6) 

where k ~ {q — 2)es(me/TOp)'~^/^e^^, with q — 2.2 as 
observed in the TW object and Cmin = Tmin'TieC^, one can 
obtain the total synchrotron luminosity 



14.7 GHz 



Ls 



Ls{v)du = 10-^/3/i7™„M_6«8i?3So.oi Lq, (7) 



where i?o.oi = B/(Q.Q1 G). The synchrotron luminosity 
of the TW object in the frequency range v = 1.6 — 14.7 
GHz is L = 5 X 10^® i©, as inferred from Reid et. al 
1)1995(1 . It is seen that at fiducial values of the parame- 
ters ~ 5 X 10~37~°„^ provides a synchrotron luminos- 
ity from the boundary layer comparable to the luminosity 
observed in the TW object. Note, that with the same as- 
sumptions about the generation of relativistic electrons in 
the jet with a coUimation angle h its total synchrotron 
luminosity is Lsj = hLso ^ 1, so that the observed 
luminosity of the TW object can be understood only if 
[3h? ~ 5 X 10~'^7~^„^, i.e. under the jet hypothesis a much 
larger (by a factor I//1 ^ 1) fraction of the wind energy 
should be transformed into the shocks which accelerate 
electrons. 



3. Angular size versus frequency 



Reid et al. H1995|l found that the angular size of the syn- 
chrotron source decreases with frequency, and attributed 
it to self-absorption of synchrotron photons. In our model 
the maximal (along radius) optical depth is 



2 / ^(r)dr. 



(8) 



where is the inner disk radius, ^ is the absorption co- 
efficient for synchrotron emission (see Ginzburg 1975) 



3e 



2T:me V 2Trm^c^ 



q/2 



^5(9+2)/2^-(9+4)/2 



(9) 
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g{q) ~ 0.7 for g ~ 2 (see Ginzburg ll975|l . This gives for the 
above assumption about the energy density in relativistic 
electrons e{r) = /?Muoo/87rr^, and for B ~ i3o(?'o/^)'' 
optical depth 



one obtains 



r = 3 X 10**^B^''o+o^/^M_67-o \ at = 10 GHz, 



and 



r = 4 X 10"/5i3(^otopM_6ro 1, 



ai V — 1 GHz, 



(10) 



(11) 



and obviously r is much less than one for ro > R*. 
One can show that free-free absorption is also negligible 
for the parameters of interest. Note that for the best fit 
parameters derived by Reid et al. within the jet model 
(cmm = 15meC^, the total density of relativistic elec- 
trons / N{e)de = 5 x lO"^ croT^ , Bq = 0.01 G) Eq. 
® gives optical depth r lO^^^L < 1 at = 1 GHz, 
L being the characteristic size of the absorbing region. 
Under such conditions the decrease of the angular size of 
the source with frequency (approximately as 9 (x ^ 
Reid et. al. 1995) can be explained only by assuming 
that the spectral properties of the synchrotron emission 
vary with radius. For a power-law spectrum of relativistic 
electrons as given by Eq. (|SJ) in the energy interval e — 
[Gmin,emax] the synchrotrou emissivity is approximately 
a power-law cx zy-(9-i)/2 in the frequency interval 
V = 0.29[vrmn, i^max], whcrc Vrmn-max corrcspoud to min- 
imal and maximal energies of electrons. Below and above 
this interval the emissivity varies as /i/ (x (j^/t'mm)^^'~^^^^, 
and oc {i^/i^max)''''~'^^^'^ exp(-iy/iy^ax), respectively. 
The mean frequency in the spectrum is determined by 
g, lymin and Vmax, and j> = (3 - q)vmaxl (5 - g) for g < 3 
and Vmax S> Vmin- The observed size of a synchrotron 
source will then decrease with frequency if the mean fre- 
quency in the spectrum v decreases with radius: at fixed 
frequency v the main contribution to the observed flux 
comes from a region with radius corresponding to the con- 
dition V < Q.29Vrnax{l')- 

For fixed boundary energies emin and emax in the spec- 
trum © the dependence of the critical frequencies Vmin 
and Vmax on radius is determined by the magnetic field 
B{r). The latter can be taken from observations of B in 
H2O maser sources by Fiebig & Glisten il989| (see also 
Liljestroem & Gwinn 2000 for more recent measurements) 



B ~ 60 



nT 



1/2 



iG. 



(12) 



For a standard a-disk (see Pringle [l981|l with temperature 
r>i?*, (13) 



T 



/ ■ \ 1/4 

/ SGMMd \ 



and kinematic viscosity 



a < 1, 



(14) 



P - 



(15) 



here M is the mass of a central star, Md, the accretion 
rate, Cg is the sound speed in the disk, i?, is the star radius, 
the other notations have their usual meaning. When T 
from Eq. H13|l and p from Eq. H15|l are substituted in Eq. 
one obtains b = 21/16 



Socr-21/16. 



(16) 



Following the above arguments we arrive at r oc j/'^^/^i 
for Craax =const. This is slightly weaker than observed in 
the TW object, however, given the observational uncer- 
tainties and the rough character of these estimates, the ob- 
tained correspondence between theoretical and observed 
size-frequency dependence can be considered as satisfac- 
tory. One should mention in this connection that a possible 
deviation of S(r) from that predicted for standard steady 
a-disks (see. Nelson, Benz, & Ruzmaikina I2()0()| . and ra- 
dial variations of the electron spectrum may contaminate 
this correspondence. It is worth noting in this regard, that 
when relativistic electrons are produced by a jet impinging 
on a dense ambient gas, an increase in surface brightness of 
the synchrotron emission is expected in the regions where 
the jet produces bow shocks and electrons are efficiently 
accelerated (in our case, these are the east and west ends 
of the synchrotron structure). When being accelerated in 
the bow shocks relativistic electrons can diffuse and fill 
the jet itself; however, their energy decreases and thus the 
surface brightness drops at lower radius (see for discussion 
Begelman, Blandford & R,ees ll983| . In such conditions, an 
inverse dependence of the angular size on frequency, i.e. 
the angular size increasing with frequency, should hold. 

4. Synchrotron structure around TW object: disk 
versus jet 

4.1. Dust emission: disk 

Wyrowski et al. (1999) found that the synchrotron source 
in the vicinity of the TW object is associated with a ther- 
mally radiating dust cloud which also has a E-W elongated 
morphology. The total luminosity of the cloud in the whole 
range of wavelengths 4 x is estimated by them 

from an upper limit on the contribution of the dust emis- 
sion of 0.5 Jy at 220 GHz and the distribution of the ro- 
tation temperature of HNCO. It is obvious that such a lu- 
minosity cannot stem from interaction of the stellar wind, 
whose mechanical luminosity is ~ 150 Af_6w| Lq, with 
the disk. It is worth noting, however, that the observed 
thermal dust emission in the 220 GHz range (~ O.SLq at 
the assumed distance to the TW object of 2.2 kpc) is only 
a small fraction ('--^ 2 x 10"'^) of the mechanical luminosity. 

The most striking finding they mentioned is that one 
of the three brightest spots of dust emission is coincident 
with the center of the synchrotron source. Wyrowski et 



4 



Shchekinov & Sobolev; Synchrotron from circumstellar disks 



al. (1999) argue that this coincidence confirms that the 
hot dust originates from interaction of the "synchrotron 
jet" suggested by Reid et. al. (1995) with the dense core 
gas confining the jet. In the disk model the observed sim- 
ilarity in morphologies of the synchrotron and dust emis- 
sion, and particularly, the coincidence of their geometric 
centers, may indicate that both are connected with ge- 
ometry of the dominant disk-like flow. The E-W align- 
ment of the dust emission and the synchrotron source, the 
three brightest dust spots (A, B and C in Wyrowski et 
al. 1999 nomenclature), and molecular emission (such as 
CH3CH2CN, CH3OH, HCOOCH3, H2CO, and SO2, see 
Fig. 3 in Wyrowski et. al. 1999) can be understood within 
a picture of a circumstellar disk embedded in a parent 
cloud and rotating with angular velocity perpendicular to 
this direction. 

4.2. Dust emission, jet 

Reid (private communication) mentions, however, another 
possible explanation: the dust giving rise to spot A (the 
brightest of the three) is heated by radiation of the mas- 
sive star producing the synchrotron jet, while the other 
two (B and C, which are located approximately on the 
continuation of the synchrotron emission to the west) are 
connected with the other massive stars. Note, however, 
that the hypothesis of the A, B and C spots arising due 
to energy radiated by embedded stars does not contradict 
the disk model itself. 

4.3. Kinematics of water masers: disk 

Fiebig et al. H1996(l have shown that the maser kinematics 
admits an ambiguous interpretation: it can be explained 
both within the disk and the outflow hypothesis. We argue 
that this ambiguity exists in the case of the TW object as 
well. 

Alcolea etal. (1993) showed that the distribution of the 
maser spots in W3(0H) consists of 3 major clusters - one 
in the center plus eastern and western clusters located at 
approximately the same distance from the center. This is 
reminiscent of the masers which appear in the disk with 
the bright source in the center and a central decrease in the 
distribution of masing water. Such a situation in the close 
vicinity of the TW object is quite plausible. However, care- 
ful analysis of the situation requires thorough modelling 
which takes into account special kinematics of the turbu- 
lent disc and changes in the maser pumping conditions. 
This is beyond the scope of the present paper. 

One possibility to tell the difference between the two 
hypotheses comes from the existence of the spectral fea- 
tures with velocities greatly different from the systemic 
velocity of the object: the maser velocities spread from -93 
km/s to -17 km/s (Cohen . 1979,1 while the systemic veloc- 
ity is about -48 km/s. Recent sensitive observations of the 
water masers in W3(H20) displayed features for which the 
difference in velocities with the bulk of the TW material 



reaches 70-90 km s~^ (Sobolev et al., private communica- 
tion). The high velocity maser features can appear in the 
clumps which are torn off the surface close to the central 
parts of the disk. In this case maser spots corresponding 
to blue-shifted and red-shifted features are likely to have 
similar distribution in the field of view. Proper motions 
of the maser features under the disk hypothesis are deter- 
mined by the proper motion of the TW object itself, the 
disk rotation and the turbulent motions at the disk sur- 
face. The procedure of estimating the TW object proper 
motion should be the same under both hypotheses and it 
is reasonable to assume that it is about 20 km/s as was 
estimated by Alcolea etal. (1993). The sense of the disk 
rotation can be guessed from the velocity gradient of the 
highly excited molecules in closest vicinity of the TW ob- 
ject observed by Wyrowski etal. (1997): it should be east- 
ward for the blue-shifted and westward for the red-shifted 
features. The value of the rotation velocity for the mate- 
rial responsible for the high velocity features is uncertain. 
However, in the central vicinity of the disk this quantity 
can be considerably greater than both the proper motion 
of the TW object and the spread of the turbulent veloci- 
ties. In this case the eastward direction of the blue-shifted 
maser spots and the westward direction of the red-shifted 
ones should be pronounced. The proper motion pattern in 
the case of the jet is discussed in the next subsection and 
is greatly different. 

4.4. Kinematics of water masers: jet 

Alcolea etal. (1993) considered water maser kinematics 
under the hypothesis of bipolar outflow and found no 
major contradictions. Wc further consider the situation 
with high velocity maser features in the case of a jet. As 
was mentioned in the previous subsection, the possibility 
to discriminate between the two hypothesis can be pro- 
vided by the interferometric measurements of positions 
and proper motions of the high velocity water maser fea- 
tures. 

The results of Alcolea etal. (1993) show agreement 
with a considerable increase of the water maser veloc- 
ities toward the center. In this case the high velocity 
features are expected to be situated close to the center. 
Observations of Wilner et al (1999) show that the syn- 
chrotron feature wiggles. Anyhow, the observed jet con- 
volution is not high and the maser kinematics shows that 
the central parts of the jet are most likely moving close 
to the tangential direction. In this case the proper motion 
of each spot producing high velocity maser feature should 
be much greater than the difference between the radial 
velocity of the feature and that of the TW object. These 
differences for the high velocity features considerably ex- 
ceed the TW proper motion (about 20 km/s according to 
Alcolea etal., 1993). So in contrast to the case of the disk 
the proper motions of blue-shifted spots which trace the 
jet should be always great and should have different direc- 
tions: westward for the spots to the west of the center and 



Shchekinov & Sobolev: Synchrotron from circumstellar disks 



5 



eastward for the spots to the east of the center. Further, in 
contrast to the case of the disk, the proper motion pattern 
for the red-shifted spots which trace the jet should be the 
same as that for the blue-shifted ones. 

4.5. Wiggling: disk 

Wilner et al (1999) have found that the synchrotron struc- 
ture is not precisely straight, but rather wiggles at the 
WE ends with a remarkable point-like; symmetry around 
the center. In the disk model a wiggling synchrotron struc- 
ture may reflect a warped structure of the disk itself, which 
can be produced either by gravitational perturbation from 
companion stars or sufficiently massive gas clumps, or by a 
large-scale hydrodynamical instability driven by radiation 
from the central star (Pringle 1996, Armitage & Pringle 
1997). 

Wyrowski et al. (1999) estimate the star embedded in 

the dust spot C to be a BO star, and a total mass of dust 
emitting gas of ~ 15Mq. The corresponding gravitational 
perturbation on the disk from both the star and the gas 
in the B and C clumps can be evaluated by comparing 
the vertical components of the gravitational force from the 
central star ~ MH/r^, where H is the vertical scale height 
of the disk, and from a companion star and surrounding 
gas ~ (Mc -|- MBc)a/r^, where a is the distance between 
the disk plane and the center of mass of clumps B and C 
and an embedded star. The geometry of the TW source 
and of dust spots A, B and C suggests that a ^ H, which 
leads to the conclusion that the perturbation from dust 
clumps B, C and the star is sufficient to warp the disk if 
(Mc + Mbc) ~ M. 

A complementary source for keeping the disk warped 
may be connected with the interaction of the light from the 
central star with the disk: when the intrinsic luminosity of 
a central star settling onto the main sequence overwhelms 
the total disk luminosity, it can initiate a large scale in- 
stability and bend the disk (Pringle 1996, Armitage & 
Pringle 1997). The threshold luminosity for the instabil- 
ity to grow is ^ IOLq, and is much less than the total 
luminosity of the stars embedded in the TW structure, 
Ltw ~ 3 X lO^ig (Wyrowski et al. 1999). The warp ex- 
cited initially at the inner disk diffuses outside and can last 
tens of Myrs. At later times, t = 15 Myr, the tilt angle be- 
tween the inner (r ~ 10 AU) and outer (r ~ 100 AU) disk 
can reach ^ 10° (Armitage & Pringle 1997). Armitage & 
Pringle argue that this mechanism is most likely respon- 
sible for warping of the inner disk in /3 Pictoris (Burrows 
et al. 1995). As shown by Quillen (2001) a strong stellar 
wind is also expected to produce disk warping. 

A change of the wiggle morphology in the disk model 
can appear either due to relative motions of the central 
star (traced presumably by the A spot) and its com- 
panions (traced by B and C spots) if warping is pro- 
duced by gravitational perturbation from these sources, 
or by proper rotation of the disk. In both cases the char- 
acteristic times are much longer than in the jet model: 



t - s/GMjr^ ~ 5 X 10^ yr for M ~ lOM©, and r ~ 0.01 
pc. This is consistent with observed lack of morphology 
variations over 510 days (Wilner et al. 1999). 

4.6. Wiggling: Jet 

Wilner et al (1999) argue that wiggling of the synchrotron 
source in the TW may be the result of precession of the 
jet axis. In principle, precession is typical for jets in the 
vicinity of binary stellar objects and in AGNs. In the TW 
object precession can be maintained also by the gravita- 
tional influence from a companion star associated with 
the B and C dust spots. Large scale wiggling can be con- 
nected with a long-periodic precession with characteristic 
time t ^ L/u, where L is the length of the jet, and u is 
the speed of flow within it. The synchrotron structure in 
the TW object has an extent 2L ~ 0.01 pc (Reid et al. 
1995), and the velocities of the water masers provide a 
lower estimate of the flow velocity: w ~ 20 km s""'^ at 1" 
and increasing to the center (Alcolea et al. 1993), u ^ 40 
km (Cohen 1979). Recent sensitive observations of the 
water masers in W3(II20) displayed features for which the 
difference in velocities with the bulk of the TW material 
reaches 70-90 km s~^ (Sobolev et al., private communi- 
cation). Thus, the upper limit for long-periodic variations 
of the wiggle morphology is 250 yr, while the lower limit 
can be not less than five years for the jet flow slower than 
1000 km s~^ expected for the conditions in the vicinity of 
a young OB star. These estimates do not contradict the 
data of Wilner et al. (1999) on the absence of noticeable 
variability within 510 days. 

4.7. Polarization: disk 

The jet and the accretion disk surface are expected to 
show different polarization properties. In the first case 
the emission is expected to be linearly polarized because 
it arises in the regions with an apparently stronger and 
rather regular magnetic field of the jet. In principle, the 
interface between the jet and the ambient medium is sub- 
ject to the Kelvin-Helmholtz instability, which forces the 
interface to settle into a regime with highly developed tur- 
bulent motions causing depolarization. Indeed, numerical 
simulations show that super- Alfvenic jets undergo KH in- 
stabilities; however, only the relatively narrow region close 
to the jet walls becomes involved in a strongly nonlin- 
ear turbulent regime, while the body of the jet remains 
more regular since the characteristic growth rate of the 
instability there is an order of magnitude smaller than at 
the interface. As a result, the modelled synchrotron emis- 
sion reveals a high degree of polarization (Hardee & Rosen 
1999). 

Wilner et al. (1999) mention that morphologically the 
synchrotron structure in the TW source is similar to that 
observed in AGNs, where polarization in the jet compo- 
nent can reach 10%. In some cases the degree of polariza- 
tion in extragalactic jets can be as high as 50 % - an up- 
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per limit for optically thin synchrotron emission (Perlman 
et al., 1999), which is consistent with the predictions of 
Hardee & Rosen (1999). 

In the case of synchrotron radiation produced in the 
boundary layer of the disk, where the acceleration of non- 
thermal electrons is intimately connected with a highly 
developed turbulence, the emission seems to be depolar- 
ized. Wilner et al (1999) reported that the synchrotron 
source in the TW docs not show polarization higher than 

0. 05 of the total intensity of the core of the synchrotron 
structure, and from this point of view it is consistent with 
the disk model. 

4.8. Polarization: Jet 

However, Wilner et al. bring forward several arguments 
explaining the nondetection of polarization from the TW 
synchrotron source. The fact that the core component does 
not show significant polarization might reflect the similar- 
ity between this source and the jets in ANGs where the 
core component emission shows typically less than a few 
percent of the total emission. Note, that in the case of 
the TW source it can be connected with a short periodic 
(a few days) variations of the core component produced 
by a strong influence of the magnetic field of the mag- 
netosphere of the central star on the foot region of the 
jet, and highly developed random motions caused by this 
(see, e.g. Lovelace, Romanova, & Bisnovatyi-Kogan 1995, 
Goodson, Bohm, & Winglee 1999). The jet component, 
where AGN jets normally show polarization, in the case 
of the TW source is too weak for the confident detection 
of polarization. 

On the other hand, internal Faraday depolarization 
produced by thermal electrons with rie ^ 10'^ cm""^ would 
be sufficient to explain the low level of polarization. The 
upper limit derived by Wilner et al. for rig from the lack 
of detected thermal emission in the cm- wavelength range 
and from the size of the TW synchrotron structure is more 
than two orders of magnitude larger than this value. Under 
these circumstances a more precise determination of the 
thermal electron density would be of critical importance. 

5. Conclusions 

1. We have shown that synchrotron emission observed 
from the TW object can, in principle, be con- 
nected with the accretion disk around a massive star. 
Interaction of the wind and the strong shocks devel- 
oped in the wind flow with the disk surface can gener- 
ate multiple shocks in the boundary layer on the disk 
surface, which in turn can efficiently produce relativis- 
tic electrons. If only as little as than 0.5% of the wind 
kinetic energy is transformed into relativistic electrons, 
the total synchrotron luminosity can be comparable to 
that observed around the TW object in the W3(0H) 
region. 

2. The observed dependence of the angular size of the 
source on frequency can be explained as due to ra- 



dial variation of the upper cutoff frequency of the syn- 
chrotron spectrum, which can be connected in partic- 
ular with the radial dependence of the magnetic field 
in the disk. 

3. We have shown that existing data on the spatial distri- 
bution of dust emission, kinematics of molecular ma- 
terial as well as morphology and polarization of the 
synchrotron emission around the TW object can be 
explained both with the disk and the jet hypothesis. 
Further observational efforts in the directions shown 
(e.g., measurement of proper motions of the water 
maser spots corresponding to the high velocity fea- 
tures) can distinguish between the two hypotheses. 

In this paper we did considered synchrotron emission 
only around the TW object. Recently there have been re- 
ports of other sources of combined radiocontinuum and 
water maser emission associated with massive young stel- 
lar objects which probably contain disks (see Menten & 
van der Tak, 2003 and references therein). The bulk of the 
radiocontinuum emission around these sources looks dif- 
ferent from synchrotron emission as can be inferred from 
their spectral indices. Thus, the suggested phenomenon 
can be widespread in the massive objects but difficult to 
observe because of confusion with other types of emission. 

In our model the synchrotron luminosity of the cir- 
cumstellar disk is proportional to the mass loss rate M, 
the wind velocity v, and the disk radius R; it depends 
on the magnetic field slightly more strongly, as B^'^. One 
can thus think that only the disks around massive stars 
with high enough M and v can be sufficiently bright in 
synchrotron radiation. Note, howc^Tr, that among the es- 
sential factors which can determine synchrotron luminos- 
ity are also the separation between the shock waves Is, 
and the spectrum of the turbulent magnetic field ^ k^^; 
additional study is needed to draw firm conclusions on 
whether disks around less massive stars, such as for in- 
stance T Tauri, can shine in synchrotron radiation. 
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